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Abstract: 


Ultra high performance concretes (UHPCs) are cementitious composite materials with high level of perfor¬ 
mance characterized by high compressive strength, high tensile strength and superior durability, reached by low 
water-to-binder ratio, optimized aggregate size distribution, thermal activation, and fiber reinforcement. In the 
past couple of decades, more and more UHPCs have been developed and found their ways into practice. Thus, 
the demand for computational models capable of describing and predicting relevant aging phenomena to assist 
design and planning is increasing. This paper presents the early age experimental characterization as well as the 
results of subsequent simulations of a typical UHPC matrix. Performed and simulated tests include unconfined 
compression, splitting (Brazilian), and three-point-bending tests. The computational framework is formulated by 
coupling a hygro-thermo-chemical (HTC) theory and a comprehensive mesoscale discrete model with formulated 
aging functions. The HTC component allows taking into account various types of curing conditions with varying 
temperature and relative humidity and predicting the level of concrete aging. The mechanical component, the 
Lattice Discrete Particle Model (LDPM), permits the simulation of the failure behavior of concrete at the length 
scale of major heterogeneities. The aging functions relate the mesoscale LDPM mechanical properties in terms 
of aging degree, defined in this work as the ratio between the quasi-static elastic modulus at a certain age 
and its asymptotic value. The obtained results provide insights in both UHPC early age mechanisms and a 
computational model for the analysis of aging UHPC structures. 

1 Introduction 

Ultra high performance concretes (UHPCs) are cementitious composites characterized by high compressive 
strength, typically greater than 120 MPa (17 ksi), low water-binder ratio, optimized gradation curve, use 
of thermal activation, fiber reinforcement and superplasticizers. Moreover, UHPC has a discontinuous pore 
structure that reduces liquid ingress and permeability [Ij, leading to significantly enhanced durability, longer 
service life and lower costs for maintenance. UHPC became commercially available in the beginning of the 21st 
century and has been utilized in the construction industry especially for bridge applications and tall buildings 
around the world across North America, Europe, and Asia. 

Among various types of UHPCs developed by researchers around the globe. Ductal is the most commonly 
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used and is commercially available in the United States, Europe and Asia. Other types of UHPCs include but 
are not limited to compact reinforced composite (CRC) developed by Aalborg Portland in 1986, the UHPC mix 
of Teichmann and Schmidt in Germany, and CorTuf developed by the U.S. Army Corps of Engineers Research 
and Development Center (ERDC), which all utilize water/cement ratios within the range of 0.15~0.21, fine 
aggregate gradation, silica fume, superplasticizers and steel fibers [2J. The inclusion of silica fume allows a 
pozzolanic reaction to produce additional calcium-silicate-hydrate (C-S-H), which can fill the would-have-been 
voids and provide a much denser matrix. The perfectly graded fine aggregate plays a similar role but as filler in 
minimizing the number and volume of possible voids. 

While more and more UHPCs are developed and utilized in the construction industry, what is lacking in the 
available literature is a model for the evolution of material properties of UHPCs at early age. This is crucial 
in terms of structural design, project planning, and building optimization. The construction of bridges and 
skyscrapers, for example, structures in which UHPC materials are utilized the most, are typically characterized 
by a complex sequence of construction stages which may last for months or even years. Being able to model 
how the material ages and at what pace each mechanical property grows could greatly help optimize the 
construction sequence and hereby reduce the total construction time and budget. The behavior at early age is 
also of relevance for the long-term performance and durability due to e.g. restraint forces [3], or cracking at 
early age |4l|5]. Consequently, such a computational framework is quintessential to the accurate performance 
assessment of degrading concrete structures since all long-term analyses are anchored to the construction stages 
[ 6 ]. 

Moreover, the demand for computational models that can be utilized in design is growing with the increased 
adoption of UHPC materials in practice. For example, Chen and Graybeal [7J investigated the behavior of 
existing UHPC structural components including prestressed girders by finite element methods and showed good 
abilities of replicating the structural response. However, conventional constitutive formulations are unable to 
capture the evolution of material properties and, correspondingly, the time-dependent structural response of 
concrete as an aging material. This is an important aspect to consider given that the current construction 
practice requires shorter and shorter construction times associated with, for example, early load applications 
and prestressing. While the contributions of fibers are quite essential for the success of UHPCs, this paper 
focuses solely on a computational framework for the evolution of the matrix properties. The presence of fiber 
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reinforcement is not directly investigated here. However, related research efforts have shown that the behavior 
of fiber reinforced UHPCs can be predicted based on a calibrated constitutive model for the plain concrete and 
an appropriate fiber concrete interaction model for the discretely placed fibers [8J. 

It is well known that the strength of cement based composites such as concrete increases rapidly at early 
age. However, the chemical and physical mechanisms behind this phenomenon are complex and consist of 
multiple coupled components. The cross-effects between hydration reaction, temperature evolution, and member 
deformation involve complex chemo-physical mechanisms that operate over a broad range of length and time 
scales, from nanometer to meter, and from fractions of seconds to years [9J. Notably, evolution laws for maturing 
concrete based on Arrhenius-type time acceleration concepts are widely supported by a good agreement with 
experimental data [TPI [TTl IT^. 

Ulm and Coussy [12] studied the thermo-chemo-mechanical coupling of concrete at early age with a for¬ 
mulation based upon thermodynamics of open porous media composed of a skeleton and several fluid phases 
saturating the porous space. It accounts explicitly for the hydration of cement by considering the thermo¬ 
dynamic imbalance between the chemical constituents in the constitutive modeling at the macrolevel of the 
material description. However, the effects from stress and temperature evolutions were neglected. Afterwards 
they extended the thermo-chemo-mechanical cross effects characterizing the autogenous shrinkage, hydration 
heat and strength growth, within the framework of chemoplasticity [9J. Cervera et. al. [T3j applied the reactive 
porous media theory and introduced a novel aging model which accounts for the effect of curing tempera¬ 
ture evolution featuring the aging degree as an internal variable. They suggested that the evolution of the 
compressive and tensile strengths and elastic moduli can be predicted in terms of the evolution of the aging 
degree Ha ttSj. The model considers the short-term mechanical behavior based on the continuum damage 
mechanics theory and the long-term mechanical behavior based upon the microprestress-solidification theory 
[16J. Bernard, Ulm and Lemarchand [I7| developed a multi scale micromechanics-hydration model to predict 
the aging elasticity of cement-based materials starting at the nano level of the C-S-H matrix. Lackner and 
Mang [18J proposed a 3-D material model for the simulation of early-age cracking of concrete based on the 
Rankine criterion formulated in the framework of multi surface chemoplasticity. Gawin, Pesavento, and Schrefler 
[I^ [20] proposed a solidification-type early-age model and extended it to account for coupled hygro-thermo- 
chemo-mechanical phenomena, which was applied to solving practical problems, such as analysis, repair, and 


4 


rehabilitation of concrete structures [211122]- 

Di Luzio and Cusatis [221124] formulated, calibrated, and validated a hygro-thermo-chemical (HTC) model 
suitable for the analysis of moisture transport and heat transfer for standard as well as high performance 
concrete. In this study, classical macroscopic mass and energy conservation laws were formulated with humidity 
and temperature as primary variables and by taking into account explicitly various chemical reactions including 
cement hydration, silica fume reaction, and silicate polymerization [23]. Furthermore, Di Luzio and Cusatis [25] . 
amalgamated the microplane model and the microprestress-solidification theory. This unified model takes into 
account all the most significant aspects of concrete behavior, such as creep, shrinkage, thermal deformation, 
and cracking starting from the initial stages of curing up to several years of age. 

While continuum mechanics and finite element solvers are broadly utilized for mechanical analysis of concrete 
structures at the macroscopic levels, the Lattice Discrete Particle Model (LDPM) [261 EZ] provides additional 
insights into failure behavior of concrete at smaller length scales. LDPM simulates concrete at the length scale 
of coarse aggregate pieces (mesoscale) and is formulated within the framework of discrete models, which enable 
capturing the salient aspects of material heterogeneity while keeping the computational cost manageable [26] . 

The HTC model and LDPM are selected as basis for the early age mechanical model formulated in this study. 
On the one hand, the HTC model can comprehensively capture the hygral and thermal evolutions and chemical 
reactions during aging. On the other hand, LDPM can provide insights into the concrete behavior at the 
mesoscale level and also simulate well the mechanical behavior of concrete structures under various loads at the 
macroscopic level. The two models are connected by a set of proposed aging functions that form the core of the 
proposed HTC-LDPM coupled early age framework, allowing prediction of the behavior of concrete members 
at early age and beyond. In support of the calibration and validation tasks, a comprehensive experimental 
investigation for a UHPC at early age was also carried out. The performed tests include unconfined cylinder 
cube compression, cylinder splitting, notched and unnotched beam three-point-bending tests, and internal 
relative humidity measurements on cylindrical specimens. The evolution of unconfined compressive strength, 
tensile splitting strength, tensile flexural strength, and the associated evolving stress-strain relationships are then 
extracted from the experiments and presented in terms of concrete aging. 
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2 Computational Framework for the Simulation of UHPC at Early Age 


The proposed hygro-thermo-chemo-mechanical early-age model for cement based concrete consists of three 
major components: the hygro-thermo-chemical (HTC) model, the Lattice Discrete Particle Model (LDPM), 
and the formulated aging functions bridging the two. 

2.1 Hygro-Thermo-Chemical (HTC) model 

The behavior of concrete at early age heavily depends on moisture content and temperature. The overall 
moisture transport can be described through Pick’s law that expresses the flux of water mass per unit time J 
as a function of the spatial gradient of the relative humidity h. Moisture mass balance requires the variation 
in time of the water mass per unit volume of concrete to be equal to the divergence of moisture flux, div J. 
The water content w is the sum of evaporable water We (capillary water, water vapor, and absorbed water) and 
non-evaporable (chemically bound) water Wn [2H1 ESI [30]. Assuming that We is a function of relative humidity 
h, degree of hydration Uc, and degree of silica fume reaction one can write Wg = We{h, Uc, as), which 
represents an age-dependent sorption/desorption isotherm. Consequently, the moisture mass balance equation 
reads [23] : 


V • (DhVh) 


dwe dh 
dh dt 


dwp, . dwe . 

O-c + ■?:—Oig + Wj 


= 0 


( 1 ) 


dac da 

where Dh is moisture permeability. The enthalpy balance is also influenced by the chemical reactions occurring 
at the early age. One can write, at least for temperatures not exceeding 100 °C 


V ■ (AtVT) -pct-^ + dccQT + dssQT = 0 (2) 

where = hydration enthalpy, = latent heat of silica-fume reaction per unit mass of reacted silica- 
fume, p is the mass density of concrete. A* is the heat conductivity, and q is the isobaric heat capacity of 
concrete. 

According to the thermodynamics based model proposed by Ulm and Coussy [121 EE| and later revised by 
Cervera et al. [I3j, the hydration kinetics can be described by postulating the existence of a Gibb’s free energy 
dependent on the external temperature T and the hydration extent Xc- For convenience the hydration degree 
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can be defined as a, 


c = Xc/xT’ where xT theoretical asymptotic value of the hydration extent in ideal 

hygrometric conditions. Analogously, the asymptotic degree of hydration, can be expressed as = xTtxT- 
With the two assumptions that (1) the thermodynamic force conjugate to the hydration extent, named 
the chemical affinity, is governed by an Arrhenius-type equation and (2) the viscosity governing the diffu¬ 
sion of water through the layer of cement hydrates is an exponential function of the hydration extent [T2] . 
Cervera et al. [I3] proposed the evolution equation for the hydration degree: otc = and 

^c(crc) = ^ci(^ + Oic){oi^ — where Ac{ac) is the normalized chemical affinity, Eac is the hy- 

dration activation energy, R is the universal gas constant, and rjc. Ad, and Ac 2 are material parameters. To 
account for the situation that the hydration process slows down and may even stop if the relative humid¬ 
ity decreases below a certain value, the equation can be rewritten as: dc = Ac.{ac)ldh{h)e ^vhere 

ldh{h) = [1 -|- (a — . The function I3h{h) is an empirical function that was first proposed for the definition 

of the equivalent hydration period by Bazant and Prasannan [33]. The parameters a and b can be calibrated 
by experimental data but values a = 5.5 and 6 = 4 can be generally adopted [331 HSj. 

The theory adopted for cement hydration can be as well utilized for silica fume (SF) reaction since the 
kinetics of a pozzolanic reaction can also be assumed to be diffusion controlled. Accordingly, the degree of 
SF reaction, as, is introduced, dg = As{as)e~^‘^’’/^'^, where = As\{^ + Q;s)(a“ — as)e~^‘’°‘A°‘?’, is 

the SF normalized affinity, Eas is the activation energy of SF reaction, and af is the asymptotic value of the 
SF reaction degree [23]. Eas/R = 9700 K can be generally assumed [S^. The rate of SF reaction degree is 
assumed not to depend on h since the pozzolanic reaction continues even at very low relative humidities [35] . 

Hydration degree and silica fume reaction degree can also be combined to define a total reaction degree [25] 
as 

^ ac{t)cQf + as{t)sQf 
afcQf + afsQf 

A detailed discussion of the calibration procedure for the HTC model can be found in Di Luzio and Cusatis [24] . 

2.2 Lattice Discrete Particle Model 

In 2011, Cusatis and coworkers [251 [2Z] developed the Lattice Discrete Particle Model (LDPM), a mesoscale 
discrete model that simulates the mechanical interaction of coarse aggregate pieces embedded in a cementitious 
matrix (mortar). The geometrical representation of the concrete mesostructure is constructed through the 
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following steps. First, the coarse aggregate pieces, assumed to have spherical shapes, are introduced into the 
concrete volume by a try-and-reject random procedure. Secondly, nodes as zero-radius aggregate pieces are 
randomly distributed over the external surfaces to facilitate the application of boundary conditions. Thirdly, 
a three-dimensional domain tessellation, based on the Delaunay tetrahedralization of the generated particle 
centers, creates a system of polyhedral cells (see Fig.[^ interacting through triangular facets and a lattice system 
composed by the line segments connecting the aggregate centers. The full description of LDPM geometry is 
reported in Cusatis et al. j2Sll2Z]- 



Figure 1: One LDPM Cell around an aggregate piece. 


In LDPM, rigid body kinematics are used to describe the deformation of the lattice particle system and the 
displacement jump, |uc], at the centroid of each facet is used to define measures of strain as 


_ n^[ucl _ l^[ucl _ m^luc] 

Cat —--; Ci — — - —; Cm — -- 


( 4 ) 


where ^ = interparticle distance; and n, 1, and m, are unit vectors defining a local system of reference attached 
to each facet. It was recently demonstrated that the strain definitions in Eq. [^correspond to the projection 
into the local system of references of the strain tensor typical of continuum mechanics [3^1 EZl I3B] . 

Next, a vectorial constitutive law governing the behavior of the material is imposed at the centroid of 
each facet. In the elastic regime, the normal and shear stresses are proportional to the corresponding strains: 
tN = E]\fe*j^ = Ejq{e]s[ — e%)-, Im = E^e]^ = Et^cm — e^); = Ej-e*]^ = Ex^ei — e°), where Ejsi = Eq, 

Et = olEq, Eq = effective normal modulus, and a = shear-normal coupling parameter; and e%, e^, e° 
are mesoscale eigenstrains that might arise from a variety of phenomena such as, but not limited to, thermal 
expansion, creep, shrinkage, and deterioration mechanisms (e.g. alkali-silica reaction). 

For stresses and strains beyond the elastic limit, the LDPM formulation considers the following nonlinear 
mesoscale phenomena [33 [40] |21]: (1) fracture and cohesion; (2) compaction and pore collapse; and (3) 





friction. 


Fracture and cohesion due to tension and tension-shear. For tensile loading (e^ > 0), the fracturing 
behavior is formulated through an effective strain, e* = \/e*^ + a{e*^ + and stress, t = + {t\j -|- t\)/a, 

which define the normal and shear stresses as = e^(t/e*); Im = «e^(t/e*); tL = ae}^{tje*). The effec¬ 
tive stress t is incrementally elastic (i = Eqc) and must satisfy the inequality 0 < t < abt{e,u) where 
o'bt = o'o(i^) exp [—iJo(i^)(e — eo(a.;))/(To(a;)], (x) = max{a;,0}, and tan(a;) = = tN^/a/tT, and 

-F e*j^. The post peak softening modulus is defined as Ho{uj) = Ht{2uj where Ht is the 
softening modulus in pure tension (ca = 7r/2) expressed as Ht = 2 Eq/ [itli— 1); it = ^E^Gtlol', t is the 
length of the tetrahedron edge; and Gt is the mesoscale fracture energy. LDPM provides a smooth transi¬ 
tion between pure tension and pure shear (ta = 0) with parabolic variation for strength given by (Jo(a;) = 
atrgt(^ — sin(a;) -F A/sin^(a;) -F 4q; 003 ^( 0 ;)/[2q; cos^(a;)], where r-gt = Ogjat is the ratio of shear strength 
to tensile strength. 

Compaction and pore collapse from compression. Normal stresses for compressive loading (e^ < 0) 
must satisfy the inequality —o'fec(eD,ey) <tN<0, where abc is a strain-dependent boundary depending on the 
volumetric strain, cy, and the deviatoric strain, = ctv — ey. The volumetric strain is computed by the volume 
variation of the Delaunay tetrahedra as ey = AV/SVq and is assumed to be the same for all facets belonging 
to a given tetrahedron. Beyond the elastic limit, —abc models pore collapse as a linear evolution of stress for 
increasing volumetric strain with stiffness He for —ey < Cd = KcoGco'- abc = o'co + (~ey — ecQ)Hc{rDy)\ 
Heirnv) = Hco/il + Kc 2 {rov — «ci)); o'co is the mesoscale compressive yield stress; r^y = e^/ey and Hci, 
Kc 2 are material parameters. Compaction and rehardening occur beyond pore collapse (—ey > ed)- In this 
case one has abc = crci{r-Dv) exp [(-ey - eci)Hc{rDv)/crdirDv)] and adirov) = o'co + (cd - eco)Hc{rDv)- 

Friction due to compression-shear. The incremental shear stresses are computed as Im = Et i^*M ^m) 
and ii = Et{Gi^ where e^ = ^dep/dtM, = ^dep/dtL, and ^ is the plastic multiplier with loading¬ 

unloading conditions 9 ?.^ < 0 and > 0. The plastic potential is defined as 99 = \/t\^ + t\ — abs{tN), where the 
nonlinear frictional law for the shear strength is assumed to be abs = o's-F(/io~/^oo)o'Aro[l—exp(tAr/o'Aro)]—/ioo^Tvl 
am is the transitional normal stress; po and p-oo are the initial and final internal friction coefficients. 

Each material property in LDPM governs part of the concrete behavior under loading. The normal elastic 
modulus, which refers to the stiffness for the normal facet behavior, E^, governs the LDPM response in the 
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elastic regime, along with the coupling parameter a. Approximately, the macro scale Young’s modulus E and 
Poisson’s ratios v can be calculated as E = Eq(2 + 3a)/(4 + a) and z/ = (1 — Q!)/(4 + a). Typical concrete 
Poisson’s ratio of about 0.18 is obtained by setting a = 0.25 m 

The tensile strength, at, and characteristic length, it, together define the softening behavior due to fracture 
in tension of LDPM facets [27]. The mesoscale fracture energy can be obtained by the relation Gt = itat/2Eo. 
Calibration of at and it is typically achieved by fitting fracture-related experimental data, e.g. the load- 
displacement curve of a notched three-point-bending test. 

The softening exponent, rit, governs the interaction between shear and tensile behavior during softening at 
the facet level. One obtains more ductile behavior in both compression and tension by increasing Ut, however 
the increase is more pronounced in compression than in tension. The shear strength, as, is the facet strength 
for pure shear and affects the macroscopic behavior in unconfined compression. Compressive yielding stress, 
(Jco; initial hardening modulus, Elea', transitional strain ratio, kco', triaxial parameters, and kc 2 ', and densified 
normal modulus, Ed, define the behavior of the facet normal component under compression and affect mostly 
the macroscopic behavior in hydrostatic and highly-confined compression. The initial internal friction, /iq; and 
transitional stress, a^o, mainly govern the mechanical response in compression and have no influence on the 
tensile behavior in LDPM. At the macroscopic level they mostly affect the compressive behavior at zero or low 
confinement. Detailed descriptions of effect and function of other LDPM mesoscale parameters can be found 
in Cusatis et al. [27] . 

Besides the compatibility and constitutive equations discussed above, the governing equations of the LDPM 
framework are completed through the equilibrium equations of each individual particle. 

LDPM has been utilized successfully to simulate concrete behavior under various loading conditions [251 EZl 
[4n [421 [43] . Furthermore, it was recently extended to properly account for fiber reinforcement [3 and it 
has the ability to simulate the ballistic behavior of ultra-high performance concrete (UHPC) |48j. In addition, 
LDPM showed success in structural scale analysis using multiscale methods [371 [49l [50l 1441 [45l [46]. 

2.3 Aging Formulation 

Experimental studies show that elastic modulus (and strength) evolution at early age depends not only on the 
degree of chemical reactions, but also, independently, on the curing temperature history [51] HH [52| . To account 
for this effect the aging degree A is introduced: 
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0 


fora < ao 


(5) 


A = 

( ) iBx — 2Axa)6i fora > ao 

I \J-max~J-ref J ^ ' 

where T^ax = 120 °C, T^e/ = 22 °C, Bx = [I + Ax{al^ - ao)]/(aoo - ao), aoo = 1, ao = total reaction 
degree at setting, nx and Ax are model parameters obtained from fitting experimental data related to aging 
mechanical properties, and a is the total degree of reaction as previously introduced. 

At constant curing temperature, experiments show that elastic modulus and strength are monotonously 
increasing functions of concrete age |521 [Ml [M]- Hence, the aging degree must also be a monotonically 
increasing function of the total reaction degree. This can be enforced by imposing the following restriction on 
the parameter Ax'. Ax < A^^^ = (2aoo — 2ao — a^ + a^)"^. Furthermore, it is worth pointing out that with 
typical values (0.1-0.3) of the hydration degree at setting, ao, for Ax < 0.25 and constant temperature, the 
function A(a) is practically coincident with a straight line. 

The aging degree can be calibrated either with elastic modulus data or with strength data relevant to 
various concrete ages and at least two different temperatures. In this study elastic modulus data are used and, 
consequently, the normal modulus, Eq, is assumed to have a linear relation with aging degree A: 

Eo = E^X ( 6 ) 

Tensile strength, at, compressive yielding stress, dc, and transitional stress, a^o, on the other hand, are 
assumed to have power-law type relations with aging degree: 


di = d^A'^*; de = drA’^- d^o = 


(7) 


where Ut, Uc, uno si'e positive constants. For simplicity and without sacrificing the quality of the simulations 
the power-law exponents can be equated, introducing a single aging exponent Ua = rit = Uc = nm- Lastly, the 
tensile characteristic length, it, is assumed to be a linear decreasing function of aging degree. This dependence 
is introduced to describe the well known brittleness increase with age: 


4 = £r(fca(l-A) + l) 


( 8 ) 
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where ka is a positive constant. All the aging functions are formulated such that the corresponding parameters 
approach their asymptotic values for A approaching the value of 1. The other LDPM mesoscale parameters, 
are assumed age-independent due to a lack of relevant experimental data on the response in compression under 
confinement. 


3 Experimental Characterization of Early Age Behavior of UHPC 


In order to calibrate and validate the proposed early age model, an experimental campaign was carried out 
to characterize the early age mechanical behavior of a UHPC material. The mixture proportions for the adopted 
mix design are reported in Table 1. The material composition consists of LaFarge Type H cement, F-50 Ottawa 
sand, Sil-co-sil 75 silica flour, Elkem ES-900W silica fume, ADVA-190 Superplasticizer and tap water. The 
maximum particle size, 0.6 mm, is limited to that of silica sand, which is a foundry grade Ottawa sand [56J. 


Table 1. Constituents and Mixing Proportions of UHPC CorTuf 


Ingredient 

Type 

Proportion 

Weight per kg 

Cement 

Lafarge Type H 

1.0000 

0.3497 

Sand 

F-50 

0.9674 

0.3383 

Silica Flour 

Sil-co-sil 75 

0.2768 

0.0968 

Silica Fumes 

Elkem ES-900W 

0.3890 

0.1360 

Superplasticizer 

ADVA-190 

0.0180 

0.0063 

Water 

Tap Water 

0.2082 

0.0728 


The baseline curing regime for the adopted UHPC [56J consisted of casting on day 0, demolding on day 
1, curing in a 100% humidity room (HR) with elevated room temperature of approximately 28 °C until day 7, 
followed by curing for 4 days in 85 °C water bath (WB), and 2 days drying in the oven at 85 °C. Preliminary 
unconfined compression tests using 50.8x101.6 mm (2x4 in) cylinders were carried out at nominal ages of 1 
day, 7 days, 11 days and 13 days. At least 3 specimens were tested for each age to obtain mean and standard 
deviation of material properties. Measured compressive strengths for the 4 different ages were 13 MPa ± 31%, 
66 MPa ± 23%, 123 MPa ± 29%, and 120 MPa ± 21%. As seen from the results, oven curing does not 
provide the expected increase in compressive strength. Rather, a very slight (-2%) decrease in the average 
compressive strength was observed. Considering the number of tested specimens (3) and the amount of scatter, 
no statistically significant change in mean value was detected. Consequently, later experiments disregarded the 
oven curing procedure. 
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To study the effects of hot water bath curing on strength gain, two curing protocols with and without hot 
water bath curing were explored. A first group of specimens was kept in the humidity room (HR) for 14 days, 
a second group, instead, was kept in the humidity room for 7 days after which it was placed in hot water bath 
(WB) at 85 °C for another 7 days. Both groups were later exposed to the same laboratory conditions (about 
22 °C and 50% RH). Unconfined compression tests using 50.8x50.8 mm (2x2 in) cylinders, three-point-bending 
(TPB) tests with half-depth notched 25.4x25.4x127 mm (1x1x5 in) beams, and tensile splitting tests using 
76.2x25.4 mm (3x1 in) disks were carried out on nominal ages of 3, 7, 14, and 28 days. Circumferential 
expansion control for compression tests and extensometer control for three-point-bending and splitting tests 
were utilized to avoid brittle failure and to obtain full post-peak information. Furthermore, a friction reduction 
layer (moly dry lubricant) was applied to the load platens for all compression tests. 



Figure 2: UHPC experiment setup: (a) cylinder compression, (b) splitting, (c) beam 3-point-bending 


3.1 Unconfined Cylinder Compression Test 

Unconfined compression tests were performed in a closed loop servo-hydraulic MTS load frame with a maximum 
capacity of 4.4 MN (1000 kips). A circumferential extensometer was utilized to control expansion of the 
specimens with the goal to obtain post-peak response (see Fig.|^). After centering the specimen, a preload of 
1-5% of the expected peak was applied before the actual test commenced. 

For compression tests, two phases of control were used: first stroke control followed by circumferential 
expansion control. Both control modes applied a loading rate of 0.001 mm/s. The mode of control was 
switched at estimated 30% of the potential peak load. Compression test results can be found in Table 2 and 
Fig. where stress is calculated as a = P/A (load/area) and strain as e = Ad/H (relative displacement 
between loading platens/specimen height). Utilized for the calculations of stress and strain for each specimen. 
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the measurements of height H and diameter D with standard deviations are also shown in the figures. 

As Table 2 shows, at 14 days, samples that were cured in hot water bath exhibited about 28% higher 
strength compared with the steam curing. However, after 14 more days in air (from day 14 to day 28) the 
strength gain was almost completely lost and the water cured samples were on average only 1% stronger than 
the others. This effect was, most likely, the result of both cooling and shrinkage cracking. 


Table 2. Unconfined Cylinder Compression Test Results 


Nominal Age [d] 

Actual Ave. Age [d 

Curing 

Strength [MPa 

No. of Specimens 

3 

2.9 

HR 

53.8 ± 21.7% 

4 

7 

7.0 

HR 

78.4 ± 10.6% 

4 

14 

14.5 

HR 

99.3 ± 15.6% 

4 

28 

29.8 

HR 

115.6 ±2.2% 

3 

14 

14.4 

HR + WB 

127.2 ± 14.6% 

4 

28 

30.0 

HR + WB 

116.7 ± 16.4% 

5 


3.2 Brazilian Test 

Indirect tensile splitting (Brazilian) tests were performed in a MTS load frame with a maximum axial capacity 
of about 1 MN (220,000 lbs). The testing procedure was developed according to ASTM C496. Loading blocks 
with dimensions of 15x38x159 mm (0.6x1.5x6.25 in) were used. Initially, localized failure in compression was 
avoided by the utilization of wooden support strips. However, their usage introduced an unwanted compliance 
at the beginning of the test. Therefore, aging tests had been ran without wooden strips. Contrary to ASTM 
C496 which calls for standard cylinder dimensions of e.g. 50.8x101.6 mm (2x4 in), cylindrical disks of 76.2 
mm (3 in) in diameter and 25.4 mm (1 in) in thickness were tested (see Fig. [^). This specimen geometry 
allows a relatively more stable transition to the post-peak regime. 

After the specimen was placed firmly in the MTS Load Frame and preloaded with approximately 1% of 
estimated peak to ensure contact between specimen and loading supports, the first phase of the actual test 
was carried out in load control with a loading rate of 0.1 kN/sec. As soon as 50% of the estimated peak was 
reached, the test was switched to displacement control to save experiment time and avoid sudden brittle failure. 
The displacement was taken by means of a transducer mounted on the front face of the samples (see Fig. [^). 
A relatively stable displacement rate was found to be 5x10^^ mm/sec for the highly brittle UHPC specimens. 
After reaching a stable post-peak response in displacement control, the opening rate was accelerated to save 
testing time while ensuring stability. At least 3 specimens were available for testing on each age. In spite of all 
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efforts and great precautions only for about half of the young specimens post-peak response could be obtained, 
and for none of the higher curing degrees. The Brazilian test results, including age of concrete at testing, 
curing conditions, nominal strength, number of specimens and dimension measurements, are shown in Table 3. 
The nominal tensile splitting strength is computed as ct„ = 2P/Tihd ([2xload]/[7rxheightxdiameter]) with the 
actual dimensions for each specimen. The obtained tensile strength values do not follow the expected aging 
trend. The most likely reason can be found in the almost futile attempt to obtain post-peak data by modifying 
the loading rate during the tests which might have influenced the observed peak values. 


Table 3. Brazilian Test Results 


Nominal 

Age 

Actual 

Ave. Age 

Curing 

Splitting Tensile 
Strength [MPa] 

Number of 
Specimens 

Diameter d [mm] 

Height h [mm] 

3 

3.0 

HR 

4.5 ± 10.4% 

3 

76.3 ± 0.3% 

23.2 ± 9.0% 

7 

7.0 

HR 

5.2 ± 9.1% 

3 

76.3 ± 0.2% 

24.5 ± 4.0% 

14 

14.9 

HR 

5.5 ± 0% 

1 

76.4 ± 0.0% 

23.3± 0.0% 

28 

29.6 

HR 

4.6 ± 17.1% 

3 

76.2 ± 0.3% 

24.5± 4.5% 

14 

14.6 

HR + WB 

11.4 ± 22.7% 

3 

76.3 ± 0.3% 

25.4 ± 10.4 % 

28 

29.2 

HR + WB 

9.4 ± 12.2% 

3 

76.4± 0.2% 

24.7 ± 3.3% 


3.3 Three-point-bending Fractnre Test 

Beam specimens with 25.4x25.4x127 mm (1x1x5 in) dimensions were cast for TPB tests (see Fig. [^) with 
notches of 50% relative depth. The nominal span (distance between bottom supports) was 101.6 mm (4 in). 
An extensometer sensor was glued to the bottom of the specimen with the notch between its two feet. After 
a pre-load of up to 5% of the expected peak was applied, the specimen was loaded in crack mouth opening 
(CMOD) control with an initial loading rate of 0.0001 mm/sec. This rate was increased in the late post-peak 
phase to save total testing time while ensuring a fully recorded softening behavior. The test results and nominal 
stress-strain curves can be found in Table 4 and Fig. 0 where the nominal flexural stress is obtained by 
equation a = 3PS/2BH‘^ ([Sxioadxtest span]/[2xspecimen widthxspecimen depth^]) and nominal strain by 
e = CMOD/specimen depth, the total fracture energy Gp is calculated as the area under the force-displacement 
curve divided by the ligament area. Geometry measurement statistics, which are used to calculate a, e, and 
Gp, are also included in the figures. 

In general, all experiments show the expected increase in strength with age. This is also observed for the 
TPB tests. However, also a strength decrease from 14 days to 28 days was recorded for all test types and 
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Table 4. Three-point-bending Test Results 


Nominal Age [d 

Actual Ave. Age [d] 

Curing 

Max. Nominal Stress [MPa] 

No. of Specimens 

3 

2.9 

HR 

1.79 ± 4.3% 

3 

7 

7.0 

HR 

2.12 ± 8.4% 

3 

14 

14.0 

HR 

2.65 ± 10.2% 

4 

28 

28.0 

HR 

2.55 ± 7.3% 

4 

14 

14.4 

HR + WB 

3.65 ± 2.5% 

4 

28 

28.0 

HR + WB 

3.10 ±7.8% 

4 


curing conditions, except compression with humidity room curing. This surprising result defies theory and, at 
the present, can only be explained by drying shrinkage damage, affecting those specimens that were stored in 
air after the hot water bath. 


3.4 Relative Humidity Measurements 

A quintessential input for any mechanical analysis utilizing the proposed aging framework are time-dependent 
spatial fields of reaction degrees. These are a “by-product" of the coupled chemo-hygro-thermal analysis using 
the HTC model. The two essential fields of the multi-physics problem are relative humidity h and temperature T. 
Consequently, in order to calibrate the HTC model for the specific UHPC of the case study, h andT measurement 
data are required on top of literature values for the reaction kinetics and diffusion related parameters. In this 
investigation RH measurements at the center of the specimens as well as at the specimen boundaries (ambient 
environment) were conducted with three curing routines: 14 days HR, 7 days HR -|- 7 days WB, and fully sealed 
(self desiccation) at room temperature. For each curing method, five 50.8x50.8 mm (2x2 in) cylinders were 
measured by RH<§iT sensors, see Fig. At the time of casting, a straw with one side closed by nylon mesh 
and housing a RH(§iT sensor, was vertically inserted into each specimen with the perforated end centered in 
the specimen. The other end of the straw protruding beyond the specimen surface was sealed with moisture- 
tight silicone sealant. The RH<§^T sensors remained for the whole measurement period in the self-desiccation 
specimens, which were sealed with plastic molds (Fig. |^), silicone sealant and plastic films. For the unsealed 
samples, during and after 100% RH humidity room curing, the RHcS^T sensors were taken out and put back 
in about once a day to avoid condensation while the operation time was minimized to preserve accurate RH 
measurements. Utilizing straws instead of embedding sensors directly in the specimen was motivated by the 
intent to 1) protect the sensors, 2) be able to reuse the sensors, and 3) easily check and replace sensors in 
case of malfunction. The temperature measurement is not reported because no appreciable hydration related 
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temperature increase was observed, owing to the small dimension of the specimens and the relatively high thermal 
conductivity of concrete. A time history of RH at the center of the specimens was obtained for each curing 
routine, as shown in Fig. For the sealed specimens, the RH level at the center dropped relatively rapidly 

and reached about 69% on day 40, which is much lower than that of normal concrete and high performance 
concrete [571 [5B]. This is due to the very low water/cement ratio and the low permeability of the UHPC (the 
water/cement ratio for this UHPC is only 0.2). The amount of free water decreases relatively fast as water is 
being consumed by cement hydration and pozzolanic reactions. Overall, the three parallel experiments of each 
curing regime show very consistent RH evolutions, as shown in Figs. a ) ~[^a). 



" (d) 


Figure 3: Relative humidity measurement: (a) specimens curing, (b) measurement setup, (c) plastic molds, 
(d) RH&T sensor 


4 LDPM Modeling of Early Age UHPC Behavior: Calibration and Validation 
4.1 HTC Model Calibration and Validation 

The actual HTC model calibration is an inverse analysis task. For that purpose the aforementioned 50.8x50.8 mm 
(2x2 in) cylinder specimens of RH measurements are simulated. The HTC input data include boundary con¬ 
ditions and material properties. The environmental boundary conditions are as follows: for HR curing: 28°C 
and 100% RH for 14 days then 50% RH afterwards; for HR-I-WB curing: 28°C and 100% RH for 7 days, then 
85°C and 100% RH for the next 7 days and afterwards 50% RH at room temperature; and for self desiccation 
curing: RH sealing and room temperature. The material properties and HTC model parameters used in the 
simulations are listed in Table 5, in which the data source is also reported. The internal RH measurements 
under sealed conditions are in particular used for the calibration of the self desiccation parameters. The moisture 
diffusion parameters are calibrated by humidity room curing RH measurements of 50 days and self desiccation 
RH measurements of 40 days, starting from casting, and validated by HR-I-WB curing RH measurements of 62 
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days of data, see Figs. ~ which include also the environmental RH. The simulated RH evolution agrees 
very well with those of the experimental investigations for the applied curing regimes. Figs. |^b) ~[^b) show 
the evolutions of cement hydration degree, silica fume (SF) reaction degree, and total reaction degree. The 
total reaction degree of FIR curing increases gradually and reaches about 0.89 on day 60, while the aging degree 
under HR+WB curing reaches its asymptote of 1.0, fairly early, at around 14 days of age (0.997 to be exact). 
The thermal activation results in a higher total reaction degree which also evolves with a higher rate. Figs. 
and present the spatial distributions of relative humidity and total reaction degree, respectively, after 3, 14, 
and 28 days of HR curing (14 days in 100% humidity room with elevated room temperature of 28 °C). The 
relative humidity exhibits relatively high spatial variability, especially as the boundary condition changes at the 
early ages. On the contrary, the total reaction degree shows low spatial variability of less than about 1%, due 
to the small dimension of the tested laboratory scale specimens. 


Table 5: HTC Parameters 


Material property name [unit] 

Symbol 

Value 

Source 

Density [kg/m^j 

P 

2400 

this paper 

Isobaric heat capacity [J/kg°C] 

Ct 

1100 

m 

Heat conductivity [W/m°C] 

Xt 

5.4 

pa 

Cement hydration enthalpy [kJ/kg] 

QT 

500 

m 

Silica fume reaction enthalpy [kJ/kg] 

Qf 

780 

m 

Hydration activation energy/R [K] 

Eac/R 

5490 

ra 

Silica fume reaction activation energy/R [K] 

Eas/R 

9700 

m 

Diffusivity activation energy/R [K] 

Ead/R 

2700 

PI 

Silica fume efficiency [-] 

Spe// 

0.9 

PI 

Polymerization activation energy/R [K] 

Eap/R 

6000 

pa 

Hydration parameter [h“^] 


2 X 10® 

pa 

Hydration parameter [-] 

^c2 

1 X 10“® 

pa 

Hydration parameter [-] 

Pc 

6.5 

m 

Silica fume reaction parameter [h~^] 

^sl 

5 X 10^^ 

m 

Silica fume reaction parameter [-] 


1 X 10“® 

m 

Silica fume reaction parameter [-] 

Vs 

9.5 

m 

Moisture diffusion parameter [mm^/h] 

Do/c 

1 X 10-^ 

this paper 

Moisture diffusion parameter [mm^/h] 

Di/c 

3.1 

this paper 

Moisture diffusion parameter [-] 

n 

3.9 

this paper 

Self desiccation parameter [-] 

9i 

1.5 

this paper 

Self desiccation parameter [-] 

Uc 

^vg 

0.2 

this paper 

Self desiccation parameter [-] 

us 

^VQ 

0.4 

this paper 
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Figure 4: HTC calibration of sealed specimens for self desiccation study: (a) RH evolution and (b) 
simulated reaction degrees 




Figure 5: HTC calibration by specimens under HR curing: (a) RH evolution and (b) simulated reaction 
degrees 
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Figure 6: HTC validation using specimens under HR + WB curing: (a) RH evolution and (b) simulated 
reaction degrees 



Figure 7: Relative humidity from HTC simulations under HR curing on (a) 3, (b) 14, and (c) 28 days of 
age 
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Figure 8: Total reaction degree from HTC simulations under HR curing on (a) 3, (b) 14, and (c) 28 days 
of age 

4.2 Calibration of LDPM - Aging Degree and Aging Fnnctions 

The first step in the calibration of the aging model is the calibration of the aging degree that, in this study and 
consistently with previous work of the last author |25|. is pursued on the basis of elastic modulus measurements. 
Identified values of the LDPM normal modulus for the case of humidity room curing (HR, T = 28°C) are 
plotted in Fig. as function of the total reaction degree. As shown in Fig. [^, the best fit of the data points 
can be obtained by setting = 75 GPa, ao = 0.278, and A\ = 0.2. The parameter nx is then calibrated 
by fitting the modulus data relevant to the water bath curing at 14 days. This gives nx = 0.1. Finally, 
the parameters Uq and ka are identified by fitting strength and fracture data. The optimum fit is obtained 
by setting Ua = 2.33 and ka = 22.2. It is worth observing that the adopted aging functions require that 
E/E^ = ^ which is the typical relationship between elastic modulus and strength as 

vastly adopted in the literature |59|- The resulting aging degree evolutions for WB and HR curing regimes 
are shown in Fig. [^. Similarly to the aforementioned total reaction degree, the aging degree with WB curing 
reaches its asymptotic value, fairly early, at around 14 days of age (0.95), while the aging degree of HR curing 
increases gradually and reaches about 0.84 at 60 days of age. The thermal activation results in a higher aging 
degree which also evolves at a higher rate. 
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Figure 9: (a) LDPM normal modulus vs. total reaction degree for humidity room curing (b) Aging degree 
evolution for water bath curing and humidity room curing 

All the LDPM simulations utilize a coarse-grained aggregate size in the range 2~4 mm. In Cusatis et. 
al. [2^ |2l], it was demonstrated that, although approximated, the coarse-grained simulations provide good 
accuracy with a significant reduction in computational time. A complete list of all mesoscale parameters used 
in LDPM for all ages investigated can be viewed in Table 6, where the second column lists the asymptotic 
values. HR stands for 14 days in humidity room curing, WB stands for 7 days humidity room -|- 7 days water 
bath curing regime, after which for both curing routines the specimens are stored at room temperature with RH 
around 50%. Due to the small specimen dimensions, the spatial variability of aging degree for all investigated 
ages is below about 1%, see the standard deviation values for aging degree shown in Table 6. Hence, the 
difference in the simulation results including and excluding spatial variability is negligible for the investigated 
case, see Fig. Thus, spatially averaged aging degree and corresponding material properties are utilized for 
the simulated compression, splitting and bending tests. Also because the cylinders, cylindrical disks and beam 
specimens used in compression, Brazilian, and bending tests have similar size, the same material properties on 
corresponding ages are used for all the specimen types. 

To ensure simulations utilizing spatially averaged material properties can represent the aging phenomena 
appreciably, the effect of spatial variability on each early age was also evaluated. Fig.fl^shows, for the mid-face of 
the 50.8x50.8 mm (2x2 in) cylinder on WB14/28 days, the spatial distribution of aging related mesoscale LDPM 
parameters including normal modulus, tensile strength and tensile characteristic length. A direct comparison 
between simulations with spatially constant (averaged) material properties and spatially variable properties as 
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functions of local aging degree are presented in Fig. 13 


For the investigated specimen geometries the differences 


are quite small with a maximum error in peak loads of 9%. 


Table 6: LDPM Parameters and Material Properties 


Curing& Age 

Asym 

HR3d 

HR7d 

HR14d 

HR28d 

WB14/28d 

Ave. Aging Degree 

1 

0.5582 

0.7061 

0.7762 

0.8103 

0.9541 

iStandard Deviation [%] 


±0.008 

± 0.007 

±0.072 

±1.11 

±0.007 

NormalModulus [MPa] 

75000 

41865 

52958 

58215 

60773 

71558 

DensificationRatio [-] 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

Alpha [-] 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

TensileStrength [-] [MPa] 

13.3 

3.4 

5.9 

7.4 

8.1 

11.9 

CompressiveYieldingStress [MPa] 

500 

129 

222 

277 

306 

448 

ShearStrengtliRatio [-] 

5.5 

5.5 

5.5 

5.5 

5.5 

5.5 

TensileCharacteristicLength [mm] 

10.6 

114 

80 

63 

55 

21 

SofteningExponent [-] 

0.28 

0.28 

0.28 

0.28 

0.28 

0.28 

InitialHardeningModulusRatio [-] 

0.36 

0.36 

0.36 

0.36 

0.36 

0.36 

TransitionalStrainRatio [-] 

4 

4 

4 

4 

4 

4 

InitialFriction [-] 

0.0335 

0.0335 

0.0335 

0.0335 

0.0335 

0.0335 

AsymptoticFriction [-] 

0 

0 

0 

0 

0 

0 

TransitionalStress [MPa] 

300 

77 

133 

166 

184 

269 

VolumetricDeviatoricCoupling [-] 

0 

0 

0 

0 

0 

0 

DeviatoricStrainThresholdRatio [-] 

1 

1 

1 

1 

1 

1 

DeviatoricDamageParameter [-] 

5 

5 

5 

5 

5 

5 

FractureEnergy [J/m^] 

12.47 

15.97 

26.28 

29.47 

30.12 

21.21 


Generally speaking, the age dependent LDPM model is calibrated first by TPB (notched) and cylinder 
compression test results and later validated by the other types of tests including circular disk Brazilian, cube 
compression, and unnotched beam bending tests, as well as size effect tests. A thorough discussion of age 
dependent size effect and the evolution of fracture characteristics can be found in a companion paper [60J. 
Experimental as well as simulated results of 25.4x25.4x127 mm (1x1x5 in) beam TPB (50% notched) and 
unconfined 50.8x50.8 mm (2x2 in) cylinder compression tests can be found in Fig.[^and Fig. respectively. 
Fig. [^a) and Fig. [^a) present specimen setup and typical failure crack pattern from LDPM simulations for 
TPB and compression tests. In the subfigures. Fig. [I0|(b-g) and Fig. [^b-g), the simulation and experimental 
results on the six investigated early ages are shown. Specimen dimension measurements, ultimate strength, 
fracture energy, elastic modulus, and number of specimens are reported in the figures as well. The experimental 
curves of compression tests are modified to eliminate machine compliance by matching the elastic modulus 
indirectly obtained from beam bending tests, of which CMOD measurements are not affected by the machine 
deformation or load platen setup. The strength gain and stress-strain curve shapes from simulations match those 
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from experiments very well. The drop in strength from 14 days to 28 days as documented in the experiment 
section, however, is not captured in the simulations. After 14 days of 100% RH curing for both curing routines, 
the specimens are kept at room temperature and room RH of roughly 50%. The sudden cooling after removing 
the specimens from the hot water bath likely causes some damage which is further amplified by drying shrinkage 
cracks. The presented HTC-LDPM early age model, on the other hand, is mainly focused on simulating early 
age strength gain and mechanical responses of concrete, thus drying shrinkage phenomena as well as the sudden 
cooling are not included at the moment. 

The mesoscale fracture energy, calculated as Gt = itcr'f/2Eo, is also listed in Table 6. It is interesting to 
note that by using the aging formulation presented in this paper, the dependence of the mesoscale fracture 
energy upon aging degree can be written as follows. 


Gt 


cxD —oo2 


A 


2n„ 


2E^ 



oc 





(9) 


Hence, Gt may or may not reach a maximum for A < 1 depending upon the values of and ka, which, as 
aforementioned, are positive constants. This means that the fracture energy may exhibit a non-monotonous 
relationship with the maturity of concrete as documented in the literature for several high performance concretes 
[Ml [521 IMl [64] and for the UHPC investigated in this study. The reason is that according to Hillerborg [65j the 
fracture energy is dependent on modulus, tensile characteristic length, and tensile strength which may evolve 
at different rates. Since fracture energy is not an independent material property, it must be concluded that the 
tensile characteristic length is the superior choice for the formulation of aging fracture laws, as demonstrated 
in this paper. 

Table 7 summarizes the mean values and standard deviations of the experimental as well as simulated 
nominal splitting tensile strengths for 76.2x25.4 mm (3x1 in) circular disks. Three specimens with different 
discrete particle configurations are used in the simulations for each age. As discussed earlier, the experimental 
strengths don’t quite follow the expected aging trend, which leads to a certain discrepancy when compared to 
the strength from aging-functions-based simulations. 


The 25.4 mm (1 in) cube compression tests (Fig. 14) and unnotched 25.4x25.4x127 mm (1x1x5 in) beam 


TPB (Fig. 15) tests are all conducted with 28 days old specimens, which were cured 7 days in the humidity 
room, 7 days in hot water bath, and then in air for the remaining 14 days. The nominal peak strengths of 
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the unnotched beam bending test simulations lie within the experimental scatter, see Fig. [I5|b). In both, 
simulations and experiments, the expected brittle failure mode is observed. Specimen setups, crack patterns, 
and failure types are presented in the subfigures Fig. [^b) &l [T^a). Note that the middle section of the 
beam, with its length roughly equal to its height, is discretized by lattice discrete particles while the rest of 
the specimen is a standard finite element mesh with elastic material. This configuration is utilized both in 


notched and unnotched TPB simulations, see Fig. [^and Fig. Since the material in the outer region stays 
elastic during the entire test, the replacement of LDPM material leads to significant savings in computational 
time without loss of accuracy. Generally good agreement between simulations and experiments is observed for 
indirect tensile strength, in particular for the modulus of rupture tests but also, with some limitations, for the 
Brazilian splitting tests. It is well known that both, Brazilian splitting tests and modulus of rupture tests, are 
indirect measures of tensile strength. Thus, the splitting tests on 76.2 mm disks, performed after 14 days of 
FIR+WB curing, can be also checked against simulated unnotched bending tests, cured the same way, for which 
good agreement with the experiments was already documented. In order to exclude bias due to size-effect the 
check has to be performed on equivalently sized beams, here 6.2x25.4x381 mm (3x1x15 in = height x width 
X length). The predicted modulus of rupture of 14 MPa compares reasonable well with the experimentally 
obtained splitting tensile strength of 11.4 MPa, further confirming the earlier concluded uncertainties in the 
experimental splitting tests at later ages. Further deviations between simulation and experiment arise from 
difficulties in simulating the mechanical boundary conditions of splitting tests, a well known problem. As of the 
cube compression test, the simulations match the experimental data perfectly, as shown in Fig. [T^a). 
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Figure 10: LDPM calibration by experimental results of three-point-bend tests (a) LDPM modeling (b) 
HR 3 days (c) HR 7 days (d) HR 14 days (e) HR 28 days (f) WB 14 days (g) WB 28 days 
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Figure 11: LDPM calibration by experimental results of compression tests (a) LDPM modeling (b) HR 3 
days (c) HR 7 days (d) HR 14 days (e) HR 28 days (f) WB 14 days (g) WB 28 days 
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Figure 12: Spatial fields of mesoscale LDPM parameters: (a) normal modulus, (b) tensile characteristic 
length, and (c) tensile strength, on 28 days of age under humidity room (HR) curing 
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Figure 13: LDPM simulations with and without spatial variability for compression tests 


Table 7: Brazilian Tensile Strength from Experiments and Simulations 
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Figure 14: LDPM validation by cube compression tests: (a) LDPM simulations and experimental results 
(b) LDPM modeling 
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Figure 15: 


LDPM validation by unnotched tliree-point-bend tests: (a) LDPM modeling setup and crack 


opening (b) experimental and simulated stress-strain curves 


5 Summary and Conclusions 

In recent years ultra high performance concretes have shown not only a significant rise in popularity but also 
practical relevance. Yet, a thorough understanding of the evolution of material properties at early age is 
still lacking in spite of its significance for structural design. In this paper a comprehensive numerical and 


29 































































































experimental investigation of the early age behavior of a typical UHPC is presented. The study is based on a 
large experimental campaign entailing uniaxial compression tests, tensile splitting tests, and three-point-bending 
tests at different ages and following different curing protocols, complemented by measurements of the evolution 
of internal humidity in sealed and unsealed samples. 

In order to shed light on the evolution of macroscopic material properties, an early age model is formulated 
within the framework of mesoscale discrete models. The coupled processes of moisture transport, heat transfer, 
cement hydration, and silica fume reaction are captured by a hygro-thermo-chemical (HTC) model yielding 
reaction degrees, which are then used to formulate the aging degree describing the local maturity of the UHPC. 
The local mesoscale material properties of the constitutive model, the Lattice Discrete Particle Model (LDPM), 
are obtained through rather simple aging functions, formulated in terms of aging degree. 

Based on the experimental characterization clear trends in the evolution of material properties, namely 
unconfined compressive strength, tensile splitting strength, flexural strength as well as elastic modulus for the 
investigated UHPC are noted. These increasing trends are distorted by some of the investigated curing protocols. 
In particular, the thermal activation in the hot water bath not only accelerates the curing but likely also has 
adverse affects on the material properties stemming from damage associated with the rapid cooling. Additional 
distortions in the trends may be attributed to shrinkage damage. 

Utilizing the proposed computational early age LDPM framework in furtherance of the comprehensive ex¬ 
perimental campaign, the following conclusions can be drawn: 

• It is possible to obtain post-peak response for UHPC materials, even when they are reaching a mature and, 
hence, highly brittle state if a closed loop control with a suitable and high accuracy on specimen deforma¬ 
tion measurement (e.g. crack opening for fracture or circumferential expansion for cylinder compression) 
is combined with a sufficiently slow loading rate. 

• In the present investigation thermal activation results in the expected acceleration of the chemical reac¬ 
tion and faster strength gain. However, the often reported improvement in material properties was not 
observed, on the contrary, possibly due to the negative effects of thermal shock in the aftermath. 

• The maturity of ultra high performance concrete can be predicted accurately if the coupled processes of 
moisture diffusion, heat transfer, and chemical reactions are captured in a suitable hygro-thermo-chemical 
framework, taking into account the environmental boundary conditions during curing. 
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• The effective aging degree A is a suitable parameter to quantify the maturity of concretes arising from 
cement hydration and silica-fume reaction, even following different temperature evolutions. 

• The previously postulated linear dependence of Young’s modulus on concrete maturity is confirmed. The 
analysis of the investigated UHPC shows a highly linear dependence of the local mesoscale normal modulus 
which is related to the macroscopic Young’s modulus on aging degree. 

• The evolution of local mesoscale strength-related parameters can be predicted well by monotonously 
increasing power-law type functions, also formulated in terms of concrete aging degree. The well estab¬ 
lished empirical square-root relation between Young’s modulus and macroscopic (compressive) strength 
quantities is approximately recovered for the mesoscale properties. 

• Fracture energy may exhibit a non-monotonous relationship with the maturity of concrete. The reason 
is that according to Hillerborg the fracture energy is dependent on modulus, tensile characteristic length, 
and tensile strength, which may evolve at different rates. Since fracture energy is not an independent 
material property it is not suitable for the formulation of fracture aging laws. 

• Contrary to fracture energy the tensile characteristic length changes monotonously in all investigated cases 
and shows a linear decreasing dependence on aging degree for the studied UHPC. 

• The proposed aging framework is capable of describing and predicting the non-monotonous relationship 
between fracture energy and concrete aging. 

Future work includes the extension to and validation by normal strength concretes, and the coupling with 
creep and shrinkage. This extended LDPM framework will give access to a more substantial inverse analysis of 
existing datasets and it will allow unrivaled insights into the life-time performance of concrete structures loaded 
at and exposed to the environment at early age. 
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